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The Bogoliubov dispersion relation for the elementary excitations of the weakly-interacting Bose 
gas is shown to hold for the case of the weakly-interacting photon gas (the "photon fluid") in a 
nonlinear Fabry-Perot cavity. The chemical potential of a photon in the 2D photon fluid does not 
vanish. The Bogoliubov relation, which is also derived by means of a linearized fluctuation analysis 
in classical nonlinear optics, implies the possibility of a new, superfluid state of light. The theory 
underlying an experiment in progress to observe sound waves in the photon fluid is described, and 
another experiment to measure the critical velocity of this superfluid is proposed. 



I. INTRODUCTION 

The quantum many-body problem, with its many, rich manifestations in condensed matter physics, has had a long 
and illustrious history. In particular, superconductivity and superfluidity were two major discoveries in this field. 
Although at present much is well understood (e.g., the BCS theory of superconductivity), the recent experimental 
discoveries of Bose-Einstein condensation in laser-cooled atoms raises new and interesting questions, such as 

whether the observed Bose-Einstein condensates are superfluids, or whether persistent currents can exist in these new 
states of matter. 

Historically speaking, in the study of the interaction of light with matter, most of the emphasis has been on 
exploring new states of matter, such as the recently observed atomic Bose-Einstein condensates. However, not as 
much attention has been focused on exploring new states of light. Of course, the invention of the laser led to the 
discovery of a new state of light, namely the coherent state, which is a very robust one. Two decades ago, squeezed 
states were discovered, but these states are not as robust as the coherent state, since they are easily degraded by 
scattering and absorption. In contrast to the laser, which involves a system far away from equilibrium, we shall 
explore here states close to the ground state of a photonic system. Hence they should be robust ones. 

Here we shall study the many-body problem by studying the interacting many-photon system (the "photon fluid" ) 
near its ground state. In this paper we shall explore some theoretical considerations which suggest the possibility of 
a new state of light, namely, the superfluid state. In particular, we shall derive the Bogoliubov dispersion relation for 
the weakly-interacting photon gas with repulsive photon-photon interactions, starting both from the microscopic (i.e., 
second-quantized) level, and also from the macroscopic (i.e., classical- field) level. Thereby we shall find an expression 
for the effective chemical potential of a photon in the photon fiuid, and shall relate the velocity of sound in the photon 
fiuid to this nonvanishing chemical potential. In this way, we lay the theoretical foundations for an experiment in 
progress to measure the sound-wave-like dispersion relation for the photon fluid. We also propose another experiment 
to measure the critical velocity of this fluid, and thus to test for the possibility of the superfluidity of the resulting 
state of the light. 

Although the interaction Hamiltonian used in this paper is equivalent to that used earlier in four-wave squeezing, 
we emphasize here the many-body, collective aspects of the problem which result from multiple photon-photon in- 
teractions. This leads to the idea of the "photon fluid." Since the microscopic and macroscopic analyses yield the 
same Bogoliubov dispersion relation for excitations of this fluid, it may be argued that there is nothing fundamentally 
new in the microscopic analysis given below which is not already contained in the macroscopic, classical nonlinear 
optical analysis. However, it is the microscopic analysis which leads to the new, heuristic viewpoint of the interacting 
photon system as a "photon fluid," a conception which could give rise to new ways of understanding and discovering 
nonlinear optical phenomena. Furthermore, the interesting question of the quantum optical state of the light inside 
the cavity resulting from multiple interactions between the photons (i.e., whether it results in a coherent, squeezed, 
Fock, or some other quantum state), cannot be addressed by classical nonlinear optical methods. Thus this paper 
represents a first attempt to formulate the new concept of a "photon fluid" starting from the microscopic viewpoint, 
and to lay the foundations for answering the question concerning the resulting quantum optical state of the light. 
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II. THE BOGOLIUBOV PROBLEM 



Here we re-examine one particular many-body problem, the one first solved by Bogoliubov j^,^- Suppose that 
one has a zero-temperature system of bosons which are interacting with each other repulsively, for example, a dilute 
system of small, bosonic hard spheres. Such a model was intended to describe superfluid helium, but in fact it did not 
work well there, since the interactions between atoms in superfluid helium were too strong for the theory for be valid. 
In order to make the problem tractable theoretically, let us assume that these interactions are weak. In the case of 
light, the interactions between the photons are in fact always weak, so that this assumption is a good one. However, 
these interactions are nonvanishing, as demonstrated by the fact that photon-photon collisions mediated by atoms 
excited near, but off, resonance have been experimentally observed |7j. We start with the Bogoliubov Hamiltonian 

H = Hfree + Hint 

Hfree = ^ e{p)alap 
p 

Hint = \Y^ y{i^)al+^a\_^apaq , (1) 

Kpq 

where the operators aj, and Op are creation and annihilation operators, respectively, for bosons with momentum p, 
which satisfy the Bose commutation relations 

[op, a\] 6pq and [a^, aq] = [a],, a\]^Q. (2) 

The first term Hfree in the Hamiltonian represents the energy of the free boson system, and the second term Hint 
represents the energy of the interactions between the bosons arising from the potential energy V{k). The interaction 
term is equivalent to the one responsible for producing squeezed states of light via four- wave mixing [p^ . It represents 
the annihilation of two particles, here photons, of momenta p and q, along with the creation of two particles with 
momenta p-\- k and q — k, in other words, a scattering process with a momentum transfer k between a pair of particles 
with initial momenta p and q, along with the assignment of an energy V{k) to this scattering process. 



III. THE FREE-PHOTON DISPERSION RELATION INSIDE A FABRY-PEROT RESONATOR 

Photons with momenta p and q also obey the above commutations relations, so that the Bogoliubov theory should 
in principle also apply to the weakly-interacting photon gas. The factor e{p) represents the energy as a function of 
the momentum (the dispersion relation) for the free, i.e., noninteracting, bosons. In the case of photons in a Fabry- 
Perot resonator, the boundary conditions of the mirrors cause the e{p) of a photon trapped inside the resonator to 
correspond to an energy-momentum relation which is identical to that of a nonrelativistic particle with an effective 
mass 1^,^ of m — huj/c^. This can be understood starting from Fig. [T]. 

For high-reflectivity mirrors, the vanishing of the electric field at the reflecting surfaces of the mirrors imposes a 
quantization condition on the allowed values of the z-component of the photon wave vector, = mr/L, where n is 
an integer, and L is the distance between the mirrors. Thus the usual frequency-wavevector relation 

cjik)=c[kl + kl + kl]'/\ (3) 

upon multiplication by h, becomes the energy-momentum relation for the photon 

Eip) = c[pI+pI = c[pI+pI + nW/L^Y'^ = c[pI+pI + m\^Y'^ , (4) 

where m = hnn/Lc is the effective mass of the photon. In the limit of small-angle (or paraxial) propagation, where 
the small transverse momentum of the photon satisfies the inequality 

Pi- = [pI + pI] <-Pz = hk, = hnn/L , (5) 

we obtain from a Taylor expansion of the relativistic relation, a nonrelativistic energy-momentum relation for the 2D 
noninteracting photons inside the Fabry-Perot resonator 

E{p^) = mc^ + p'i/2m , (6) 
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where m = hmr/Lc = hcu/c^ is the effective mass of the confined photons. It is convenient to redefine the zero of 
energy, so that only the effective kinetic energy, 

€{pi_)^pl/2m, (7) 

remains. To establish the connection with the Bogoliubov Hamiltonian, we identify the two-dimensional momentum 
p± as the momentum p that appears in this Hamiltonian, and the above e(pj_) as the e(p) that appears in Eq. (|l|). 



IV. THE BOGOLIUBOV DISPERSION RELATION FOR THE PHOTON FLUID 

Now we know that in an ideal Bose gas at absolute zero temperature, there exists a Bose condensate consisting 
of a macroscopic number Nq of particles occupying the zero-momentum state. This feature should survive in the 
case of the weakly-interacting Bose gas, since as the interaction vanishes, one should recover the Bose condensate 
state. Hence following Bogoliubov, we shall assume here that even in the presence of interactions, iVo will remain 
a macroscopic number in the photon fluid This macroscopic number will be determined by the intensity of the 
incident laser beam which excites the Fabry-Perot cavity system, and turns out to be a very large number compared 
to unity (see below). For the ground state wave function 5'o(A^o) with A^'o particles in the Bose condensate in the 
p = state, the zero-momentum operators uq and operating on the ground state obey the relations 

ao|*o(A^o)) = y/N^l'^oiNo-l)) 

al |*o(A^o)) = v/iVoTT|*o(A^o + 1)) • (8) 



Since Nq ^ I, we shall neglect the difference between the factors ^/No+T and vTVq. Thus one can replace all 
occurrences of ap and aj by the c-number ^/Nq, so that to a good approximation [qq, aj] « 0. However, the number 
of particles in the system is then no longer exactly conserved, as can be seen by examination of the term in the 
Hamiltonian 

Y V{K)aial,aoao « iVo ^ V{i^)aial^ , (9) 

which represents the creation of a pair of particles, i.e., photons, with momenta k and —k out of nothing. 

However, whenever the system is open one, i.e., whenever it is connected to an external reservoir of particles which 
allows the total particle number number to fluctuate around some constant average value, then the total number of 
particles need only be conserved on the average. Formally, one standard way to compensate for the lack of exact 
particle number conservation is to use the Lagrange multiplier method and subtract a chemical potential term fiNop 
from the Hamiltonian (just as in statistical mechanics when one goes from the canonical ensemble to the grand 
canonical ensemble) 

H^H' = H- y.Nop, (10) 

where Nop = X^p'^p^p total number operator, and represents the chemical potential, i.e., the average energy 

for adding a particle to the open system described by H. In the present context, we are considering the case of a Fabry- 
Perot cavity with low, but finite, transmissivity mirrors which allow photons to enter and leave the cavity, due to an 
input light beam coming in from the left and an output beam leaving from the right. This permits a realistic physical 
implementation of the external reservoir, since the Fabry-Perot cavity allows the total particle number inside the 
cavity to fluctuate due to particle exchange with the beams outside the cavity. However, the photons remain trapped 
inside the cavity long enough so that a thermalized condition is achieved after many photon-photon interactions (i.e., 
after many collisions), thus allowing the formation of a photon fluid. 

It will be useful to separate out the zero-momentum components of the interaction Hamiltonian, since it will turn out 
that there is a macroscopic occupation of the zero-momentum state due to Bose condensation. The prime on the sums 
Spi SpK' Slpq following equation denotes sums over momenta explicitly excluding the zero-momentum 

state, i.e., all the running indices p, k, q,p + k,<7 — n which are not explicitly set equal to zero, are nonzero: 

Hint = iv^(0)ajajaoao + V'(O)^ atflpajao -I- 

p 

Y ( ^(p)«p"o"pao + ^ \v{p)ala^_paoaQ + V{p)alalapa- 
p ^ 
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Here we have also assumed that V{p) = V{—p). By thus separating out the zero-momentum state from the sums in 
the Hamiltonian, and replacing all occurrences of aq and aj by \/Nq, we find that the Hamiltonian H' decomposes 
into three parts 



where 



H' = H0 + H1+H2, (12) 



Ho = lviO)alalaoao ~ lv{0)N^ , (13) 
Hi « ^'e'(p)atap + iiVo^V(p) (al^al + a^pGp^ , (14) 



where 



e'{p)^e{p) + N^V{())^N^V(p)- iJi (16) 
is a modified photon energy, and where A^o ^^nd /i are given by 

No+ < ^o\Y^alap\^o N (17) 



and 

Here Eq = {'^q\H\'^q) is the ground state energy of H. In the approximation that there is little depletion of the Bose 
condensate due to interactions (i.e., N Nq '3> 1), the first term of Eq. (pi) (i.e., Hq in Eq. (p^) dominates, so that 



and therefore that 



Eo « In^V{0) « iiW(O), (19) 



II « W(0) « NoV{0). (20) 

This implies that the effective chemical potential of a photon, i.e., the energy for adding a photon to the photon 
fluid, is given by the number of photons in the Bose condensate times the repulsive pairwise interaction energy 
between photons with zero relative momentum. It should be remarked that the fact that the chemical potential is 
nonvanishing here makes the thermodynamics of this two-dimensional photon system quite different from the usual 
three-dimensional, Planck blackbody photon system . In the same approximation, Eq. ( p^ ) becomes 

e'ip)^eip)+NQV{p). (21) 

This is the single-particle photon energy in the Hartree approximation. 

In the same approximation, it is also assumed that \Hi \ ^ \H2\, i.e., that the interactions between the bosons are 
sufficiently weak, again so as not to deplete the Bose condensate significantly. In the case of the weakly-interacting 
photon gas inside the Fabry-Perot resonator, since the interactions between the photons are indeed weak, this as- 
sumption is a good one. 

Following Bogoliubov, we now introduce the following canonical transformation in order to diagonalize the quadratic- 



form Hamiltonian Hi in Eq. (14) 



al = u^al + Vi^a.f, . (22) 
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Here m„ and are two real c-numbers which must satisfy the condition 



2 2 1 



(23) 



in order to insure that the Bose commutation relations are preserved for the new creation and annihilation operators 
for certain quasi-particles, aj^ and a^, i.e., that 



[Q;„,a^,] = (5k, and [af^a^,'] = 
We seek a diagonal form of Hi given by 



0. 



constant 



(24) 



(25) 



where a)(K) represents the energy of a quasi-particle of momentum k. Substituting the new creation and annihilation 
operators a]l and given by Eq. (^2|) into Eq. (^5|) , and comparing with the original form of the Hamiltonian Hi in 
Eq. (|lj), we arrive at the following necessary conditions for diagonalization: 



w(k)u«z;k = -NoV{k) 



f'(«)/u(«)| 

+ «'WM«)1 



Squaring Eq. ( p6| ) and substituting from Eqs. ( |27| ) and (p8[), we obtain 



(26) 
(27) 
(28) 

(29) 



where in the last step we have used Eq. (pi|). 

Thus the final result is that the Hamiltonian Hi in Eq. ( p5| ) describes a collection of noninteracting simple harmonic 
oscillators, i.e., quasi-particles, or elementary excitations of the photon fluid from its ground state. The energy- 
momentum relation of these quasi-particles is obtained from Eq. ( [29| ) upon substitution of e(K) = /2m from Eq. (|^) 



''k^NoV{k) 



(D(k) 



1/2 



(30) 



which we shall call the "Bogoliubov dispersion relation." This dispersion relation is plotted in Fig. g, in the special 
case that V{k) — V{0) = constant. (Note that Landau's roton minimum could in principle also be incorporated into 
this theory by a suitable choice of the functional form of V{k).) 

For small values of k this dispersion relation is linear in k. This feature, together with the fact that the operator 
aJ.aK in Eq. (25) describes a density fluctuation in the fluid, indicates that the nature of the elementary excitations 
here is that of phonons, which in the classical limit of large phonon number leads to sound-like waves propagating 
inside the photon fluid at the sound speed 



nm 



m J 



(-) 



1/2 



At a transition momentum Kc given by 



= 2(to7Vo1/(kc))'/' 



(31) 



(32) 



(i.e., when the two terms of Eq. (30) are equal), the linear relation between energy and momentum turns into a 
quadratic one, indicating that the quasi-particles at large momenta behave essentially like nonrelativistic free particles 
with an energy of /2m. The reciprocal of Kc defines a characteristic length scale 



Ac = 2'kTi/ Kc = Trh/mVs 



(33) 



which characterizes the distance scale over which collective effects arising from the pairwise interaction between the 
photons become important. 

Thus in the above analysis, we have shown that all the approximations involved in the Bogoliubov theory should be 
valid ones for the case of the 2D photon fluid inside a nonlinear Fabry-Perot cavity. Hence the Bogoliubov dispersion 
relation should indeed apply to this fluid; in particular, there should exist sound-like modes of propagation in the 
photon fluid. 
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V. CLASSICAL PICTURE OF SOUND WAVES IN A NONLINEAR OPTICAL FLUID 



A classical nonlinear optical treatment of a Fabry-Perot cavity which is filled with a medium with a self-defocussing 
Kerr nonlinearity (see Fig. also indicates the existence of modes of sound-like wave propagation in the nonlinearly 
interacting light. Such a nonlinear medium could consist of an alkali atomic vapor excited by a laser detuned to the 
red side of resonance. In fact, it turns out that fluctuations in the light intensity in this medium propagate with a 
dispersion relation which is identical to that given above in Eq. ( ^o|) for the weakly-interacting Bose gas. 

To derive this dispersion relation classically, we begin by considering the planar Fabry-Perot cavity shown in Fig. |^. 
Two parallel planar mirrors of reflectivity R and transmissivity T (with R + T = 1, i.e., with no dissipation) are 
normal to the z-axis and separated by a distance L. A laser beam travelling in the +z direction is incident on the 
cavity, and there results five interacting light beams in the problem. The region between the mirrors (inside the 
cavity) contains a nonlinear polarizable medium. The classical electric field obeys Maxwell's equations, written in 
wave-equation form in CGS units as 

d^E 2 1 ^^E _ 47r d^P 

5;? + VxE-^^-^^, (34) 

where E is the (real) electric field amplitude, P is the polarization introduced in the medium, and is the Laplacian 
in the transverse coordinates x and y. This equation is supplemented by boundary conditions at the two mirrors. 
Equation ( p^ ) simplifies considerably when the following assumptions arc made: 

1. The slowly-varying envelope approximation is justified, in which case we recast Eq. (^J) in terms of the field 
envelope 8. 

2. The frequency spacing between adjacent longitudinal cavity modes is much greater than 

(a) the incident laser linewidth, and 

(b) the nonlinearity bandwidth, 

allowing us to neglect the z-dependence of the field envelope (this is sometimes called the uniform field approx- 
imation). 

3. The atomic response time is much shorter than the cavity lifetime, allowing us to adiabatically eliminate the 
atomic response (i.e., the nonlinearity is fast). 

Under these reasonable assumptions the cavity's internal field envelope is governed by the Lugiato-Lefever equation 
p3[ , written here as 

at " 2fc^^^ ^ ic^«2|f + i(Ac^)f - T{£ - £d) , (35) 

where £{x,y,t) is the internal cavity field envelope amplitude, k is the longitudinal wavenumber, lo is the laser 
angular frequency, n2 is the nonlinear index inside the cavity (n « 1 -I- n2|£p), Auj — lo — oJcav is the detuning of 
the driving laser from linear cavity resonance, F = cT/2L is the cavity decay rate, and £d{x,y) is a driving laser 
amplitude. In other contexts, Eq. (|3^) is called the Nonlinear Schrodinger (NLS) equation, or the Ginzburg-Landau 
equation, or the Gross-Pitaevskii equation. The latter two of these were introduced as descriptions of superfluid and 
of Bose-Einstein-condensed systems, with a complex order parameter 5*, which here is identified with £. 
Equation ( p5| ) has the nonlinear plane-wave solution 

£ = £q exp[i{ujn2£Q + Auj)t] (36) 

when F is negligible p^ , in which case £q can be assumed real without loss of generality. Linearizing around this 
solution by substituting the form 

£ = {£o + a{x, y, t)) exp[iiujn2£o + Aw)t], (37) 
we get the following linear equation for the fluctuation amplitude (we have assumed that \a{x,y,t) \ <^ £o): 

^ ^ ^V±a + iujn2£oia + a*) . (38) 
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Here we look for a cylindrically symmetric solution appropriate for the experimental geometry (see Fig. Substi- 
tuting the trial solution 



a{p, t) — aJo{Kp)e 



PMKp)e 



-lilt 



(39) 



where Jq{Kp) is the zero-order Bessel function, p = (x^ -I- y^)^/^ is the transverse radial distance from the origin 
of a fluctuation, and K is the wavenumber of the fluctuation, we obtain the following dispersion relation for small- 
amplitude intensity fluctuations in the light filling the cavity |T5|: 



1^ 



1/2 



(40) 



where and K are the angular frequency and wavenumber respectively of the transverse sound-like mode. 

For transverse wavelengths much longer than A,, = A/ (An)^^^, where A is the optical wavelength and An = |n2| £q 
is the nonlinear index shift induced by the background beam, the transverse mode propagates with the constant phase 
velocity 



^An = 



"2 1^0^ 



(41) 



which we identify as a sound-wave velocity. This velocity is identical to the one found earlier in Eq. (^) for the 
velocity of phonons in the photon fluid, provided that one identifies the energy density of the light inside the cavity 
with the number of photons in the Bose condensate as follows: 



(42) 



where Vcav^ the cavity volume, is also the quantization volume for the electromagnetic field, and provided that one 
makes use of the known proportionality between n2 and V{Q) [ p^|Jl^ . 

In fact, the entire dispersion relation, Eq. (|o|), found above classically for sound- like waves associated with fluctu- 
ations in the light intensity inside a resonator filled with a self-defocusing Kerr medium, is formally identical to the 
Bogoliubov dispersion relation, Eq. (|3^) , obtained quantum mechanically for the elementary excitations of the photon 
fluid, in the approximation V{k) — V{Q) — constant. This is a valid approximation, since the pairwise interaction 
potential between two photons is given by a transverse 2D pairwise spatial Dirac delta function, whose strength is 
proportional to n2 jl^JT^. It should be kept in mind that the phenomena of self-focusing and self-defocusing in 
nonlinear optics can be viewed as arising from pairwise interactions between photons when the light propagation is 
paraxial and the Kerr nonlinearity is fast JT6| , pT| . Since in a quantum description the light inside the resonator is 
composed of photons, and since these photons as the constituent particles are weakly interacting repulsively with each 
other through the self-defocusing Kerr nonlinearity to form a photon fluid, this formal identification is a natural one. 



VI. AN EXPERIMENT IN PROGRESS 



We are in the process of investigating experimentally the existence of the sound-like propagating photon density 
waves predicted above for a planar Fabry-Perot cavity containing a self-defocusing (n2 < 0) nonlinear medium (see 
Fig.|). 

The sound-like mode is most simply observed by applying two incident optical fields to the nonlinear cavity: a broad 
plane wave resonant with the cavity to form the nonlinear background fiuid on top of which the sound-like mode can 
propagate, and a weaker amplitude-modulated beam which is modulated at the sound wave frequency in the radio 
range by an electro-optic modulator, and injected by means of an optical fiber tip at a single point on the entrance 
face of the Fabry-Perot. The resulting weak time- varying perturbations in the background light induce transversely 
propagating waves in the photon fiuid, which propagate away from the point of injection like ripples on a pond. This 
sound-like mode can be phase-sensitively detected by another fiber tip placed at the exit face of the Fabry-Perot some 
transverse distance away from the injection point, and its sound-like wavelength can be measured by scanning this 
fiber tip transversely across the exit face. 

The experiment employs a cavity length L of 2 cm and mirrors with intensity refiectivities of i? = 0.997, for a 
cavity finesse of roughly 1000. The optical nonlinearity is provided by rubidium vapor at 80° C, corresponding to 
a number density of 10^^ rubidium atoms per cubic centimeter. We use a circularly-polarized laser beam, detuned 
by around 600 MHz to the red side of the ^^Rb, F = 2 F' = 3 transition of the D2 line; using this closed 
transition eliminates optical pumping into the F — 1 ground state. This 600 MHz detuning of the laser from the 
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atomic resonance is considerably larger than the Doppler width of 340 MHz, and the residual absorption arising from 
the tails of the nearby resonance line gives rise to a loss which is less than or comparable to the loss arising from 
the mirror transmissions. This extra absorption loss contributes to a slightly larger effective cavity loss coefficient F, 
but does not otherwise alter the qualitative behavior of the Bogoliubov dispersion relation, nor any of the other main 
conclusions of this paper. The above criteria (1-3) for the validity of Eq. (^5|), as well as those for the validity of the 
microscopic Bogoliubov theory, should be well satisfied by these experimental parameters. An intracavity intensity of 
40 W/cm^ results in Art — 2x 10~®, for a sound speed Vs — 4.2 x 10''cm/s and transition wavelength Ac « 1 mm. For 
this intensity, Afg ~ 8 x 10^^, so that the condition for the validity of the Bogoliubov theory, A^o ^ 1, is well satisfied. 

VII. DISCUSSION AND FUTURE DIRECTIONS 

We suggest here that the Bogoliubov form of dispersion relation, Eq. ( ^o|) or (^), implies that the photon fluid 
formed by repulsive photon-photon interactions in the nonlinear cavity is actually a photon superfluid. This means 
that a superfluid state of light might actually exist. Although the exact definition of superfluidity is presently still 
under discussion, especially in light of the question whether the recently discovered atomic Bose-Einstein condensates 
are superfluids or not Q , one indication of the existence of a photon superfluid would be that there exists a critical 
transition from a dissipationless state of superflow, i.e., a laminar flow of the photon fluid below a certain critical 
velocity past an obstacle, into a turbulent state of flow, accompanied by energy dissipation associated with the 
shedding of von-Karman-like quantized vortices past this obstacle, above this critical velocity. (It is the generation 
of quantized vortices above this critical velocity which distinguishes the onset of superfluid turbulence from the onset 
of normal hydrodynamic turbulence.) 

The Bogoliubov dispersion relation (plotted earlier in Fig. ^) consists of two regimes: (1) a linear regime, in which 
there is a linear relationship between the energy of the elementary excitation and its momentum near the origin (i.e., 
for low energy excitations) corresponding to the sound-like waves, or more precisely, to the phonons in the photon 
fluid, produced by the collective oscillations of this fluid, in which the photons are coupled to each other by the 
mutually repulsive interactions between them, and (2) a quadratic regime, in which there is a quadratic relation 
for sufficiently large transverse momenta corresponding to the diffraction of the component photons, which would 
dominate when the pairwise interactions between the photons can be neglected. A crude one-dimensional model can 
give rise to an understanding of the origin of the sound-like waves in the photon fluid: Consider a system consisting of 
identical steel balls placed on a frictionless track. This system of balls is initially motionless. Now set a ball at the one 
end of the track into motion so that it collides with its nearest neighbor. The momentum transfer between adjacent 
hard spheres on this track, as they collide with one another, sets up a moving pattern of density fluctuations among 
the balls, which propagates like a sound wave from one end of the track towards the other end. Such a sound-like 
wave carries energy and momentum with it as it propagates. 

It may be asked why the classical nonlinear optical calculation gives the same result as the quantum many-body 
calculation. One answer is that one expects classical sound waves to have the same dispersion relation as phonons 
in a quantum many-body system: there exists a classical, correspondence-principle limit of the quantum many-body 
problem, in which the collective excitations (i.e., their dispersion relation) do not change their form in the classical 
limit of large phonon number. 

The physical meaning of this dispersion relation is that the lowest energy excitations of the system consist of 
quantized sound waves or phonon excitations in a superfluid, whose maximum critical velocity is then given by the 
sound wave velocity. By inspection of this dispersion relation, a single quantum of any elementary excitation cannot 
exist with a velocity below that of the sound wave. Hence no excitation of the superfluid at zero temperature is 
possible at all for any object moving with a velocity slower than that of the sound wave velocity, according to an 
argument by Landau [|l8| . Hence the flow of the superfluid must be dissipationless below this critical velocity. Above 
a certain critical velocity, dissipation due to vortex shedding is expected from computer simulations based on the 
Gross-Pitaevskii (or Ginzburg-Landau or nonlinear Schrodinger) equation which should give an accurate description 
of this system at the macroscopic level . 

We propose a follow-up experiment to demonstrate that the sound wave velocity, typically a few thousandths of the 
vacuum speed of light, is indeed a maximum critical velocity of a fluid, i.e., that this photon fluid exhibits persistent 
currents in accordance with the Landau argument based on the Bogoliubov dispersion relation. Suppose we shine 
light at some nonvanishing incidence angle on a Fabry-Perot resonator (i.e., exciting it on some off-axis mode). This 
light produces a uniform flow fleld of the photon fluid, which flows inside the resonator in some transverse direction 
and at a speed determined by the incidence angle. A cylindrical obstacle placed inside the resonator will induce a 
laminar flow of the superfluid around the cylinder, as long as the flow velocity remains below a certain critical velocity. 
However, above this critical velocity a turbulent flow will be induced, with the formation of a von-Karman vortex 
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street associated with quantized vortices shed from the boundary of the cyUnder [|19|]. The typical vortex core size 
is given by the Ught wavelength divided by the square root of the nonlinear index change. Typically the vortex core 
size should therefore be around a few hundred microns, so that this nonlinear optical phenomenon should be readily 
observable. 
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FIG. 1. A planar Fabry-Perot imposes boundary conditions which quantize the allowed values of fe^, where z is the axis 
normal to the mirrors, in units oi L, where L is the separation of the mirrors. For a plane- wave mode which propagates at a 
small angle with respect to the z axis, there arises an effective nonrelativistic energy-momentum relation for an noninteracting, 
trapped 2D photon, whose effective mass is m = Tiuj/c? (see text). 




Momentum 

FIG. 2. The energy versus momentum of an elementary excitation in the weakly-interacting Bose gas, here in the present 
case, the photon fluid. The solid line represents the Bogoliubov dispersion relation given by Eq. (^), for the special case that 
V{hi) — V{0) = constant, and the dashed line represents a quadratic dispersion relation for a noninteracting, diffracting photon 
inside the Fabry-Perot resonator. 
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FIG. 3. Fields and coordinate system in the Fabry-Perot cavity. The applied field Einc arises from a laser beam incident from 
the left. An atomic vapor excited to the red side of resonance by the incident light fills the space (the gray area) between the two 
mirrors. The presence of these atoms leads to a self-defocussing Kerr nonlinearity (corresponding to repulsive photon-photon 
interactions) inside the cavity. 




FIG. 4. Schematic of an experiment to observe the sound-like waves in a photon fluid which fills a nonlinear Fabry-Perot 
resonator. The nonlinear medium (denoted by the gray area) is an alkali atomic vapor excited by a broad laser beam (denoted 
by the broad incoming arrow) tuned to the red side of resonance. The goal is the verify the Bogoliubov dispersion relation, 
Eq. (^o|) or (^^. An electro-optic modulator (EOM) modulates the intensity of light at a radio frequency in the MHz range, 
which is then injected by means of an optical fiber tip at a single point on the entrance face of the Fabry-Perot resonator. The 
wavelength of the resulting sound-like waves can be measured by scanning in the transverse direction the tip of another optical 
fiber across the output face of the Fabry-Perot, and by measuring the phase of the modulated pick-up signal relative to that of 
the EOM modulation signal. 
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